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The conditions which determine whether a material behaves in a brittle or ductile fashion on 
mechanical loading are still elusive and comprise a topic of active research among materials physicists 
and engineers. In this study, we present results of in silico mechanical deformation experiments from 
two very different model solids in two and three dimensions. The first consists of particles interacting 
with isotropic potentials and the other has strongly direction dependent interactions. We show that 
in both cases, the excess vibrational density of states is the fundamental quantity which characterises 
the ductility of the material. Our results can be checked using careful experiments on colloidal solids. 


I. INTRODUCTION: 

Understanding the mechanical behaviour of glassy ma¬ 
terials [1-10] has engaged the attention of materials sci¬ 
entists because of both its technological ramifications and 
scientific interest [11-15]. Glasses, in some ways repre¬ 
sent an extreme limit of a supercooled liquid whose vis¬ 
cosity has increased dramatically eg. by almost 14 orders 
of magnitude with decreasing temperature within a nar¬ 
row range, 100/T) without any major change in struc¬ 
ture reaching a value of ~ lO^^Poise [16, 17]. While the 
question of why viscosity rises so rapidly has remained 
elusive [9, 18-21], in this paper, we investigate the conse¬ 
quence of this ultra slow relaxation on the failure prop¬ 
erties of amorphous materials under external load. 

In recent experiments it was found that glassy materi¬ 
als can have much larger yield stress compared to their 
crystalline counterparts with similar compositions, but 
they fail catastrophically under external loading[14, 15]. 
Thus one of the major questions which needs to be ad¬ 
dressed in order to be able to use glassy materials for 
industrial applications, is how much plastic deformation 
can the material withstand before it fails or in other 
words will the material experience brittle fracture or duc¬ 
tile fracture [22]. In Fig.I, typical failure events are 
shown for a brittle ( top panel ) and a ductile material ( 
bottom panel). One can clearly see that the morphology 
of the fracture is very different for the two cases. For the 
brittle case one has somewhat rough fracture surface, on 
the other hand for the ductile fracture, the system forms 
a neck before breaking completely. It is believed that 
for brittle materials cavities form during tension and the 
different cavities eventually percolate to form the brittle 
crack, on the other hand for the ductile materials the 
failure mechanism does not involve cavities [23]. It is im¬ 
portant to mention that the microscopic details of this 
cavitation is still not clearly understood nor the necking 
behaviour in ductile materials. 
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FIG. 1. Snapshot picture of a brittle ( top panel ) and a 
ductile (bottom panel) materials at the point of fracture. Note 
the drastic difference between the fracture morphologies for 
these two cases. 


Recent molecular dynamics (MD) simulations have in¬ 
vestigated the intrinsic relation between the properties 
of the inter-particle potential in order to understand the 
microscopic origin of brittle to ductile transition in amor¬ 
phous solids [13]. It was shown that the degree of duc¬ 
tility is intrinsically related to the growing importance 
of the contribution of plastic modes to the vibrational 
density of states of the material. Similarly Mizuno et 
al.[24] and Goodrich et al. [25] have also observed that 
the mechanical and thermal properties of a solid are di¬ 
rectly related to the density of normal modes of vibra¬ 
tion. In [25], it is found that the physics of crystalline 
solids with defects are better described by the physics of 
amorphous solids near jamming transition if the concen- 
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tration of defects are slightly large but not large enough 
to completely destroy the bulk crystalline order in the 
system. It was shown that vibrational density of states 
for these crystalline solids with defects are very similar 
to the amorphous solids near jamming transition. Their 
main conclusion is that jamming point of an amorphous 
solids is probably the better reference point than an ideal 
crystalline state to understand the physics of crystalline 
solids with defects beyond a small critical defects concen¬ 
tration. 

The aim of this paper is to investigate the connec¬ 
tion between the mechanical and vibrational properties 
of both crystalline and amorphous solids in an unihed 
framework and try to come up with a measure which can 
characterises these properties of the system irrespective 
whether the material is crystalline or amorphous in na¬ 
ture. This is performed via extensive Molecular Dynam¬ 
ics Simulation by studying different systems in which the 
mechanical behavior of the material is varied continu¬ 
ously by tuning a control parameter like the interaction 
range of the inter-atomic potential, size disorder of dif¬ 
ferent particles in the system and the quench dynamics 
of a colloidal system. 

The structure of the paper is as follows: In Sec. II we 
will describe the model systems studied, the details of 
the parameters of the model and the how the simulations 
are done for all these model systems. The numerical ex¬ 
periments and their results are described in Sec. Ill and 
finally, the main conclusions and future directions will be 
discussed in Sec. IV. 


II. MODELS AND SIMULATIONS DETAILS: 


To show that our results are generic, we report studies 
on two very different model solids. The first one consists 
of particles interacting with isotropic potentials and is a 
more or less accurate representation of a metallic glass. 
We study this model in both two and three dimensions. 
The bulk of the results presented in this manuscript re¬ 
lates to this model. We also present results for a second 
model solid which is composed of particles with strongly 
direction dependent interactions. This solid can be used 
to mimic gels, silicates and other such anisotropic ma¬ 
terials. We describe each of these models in somewhat 
more detail below. 


A. Model A: 


The model studied is a generic glass former in both 
two and three dimensions. It is a binary mixture whose 
amount of bi-dispersity of small and large particles was 
chosen to avoid any crystallization. The particles interact 


by an inter-particle potential given by 
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where r^m is the length where the potential attains its 
minimum, and Tco is the cut-off length for which the po¬ 
tential vanishes. The coefficients a, b, and C 2 i are chosen 
such that the repulsive and attractive parts of the poten¬ 
tial are continuous with two continuous derivatives at the 
potential minimum and the potential goes to zero con¬ 
tinuously at Tco with two continuous derivatives as well. 
e is the unit of energy, and /cb = 1. 


We have performed simulations where the cut off dis¬ 
tance Tco was varied from 1.2 to 2.2 keeping all other 
parameters of the interaction potential unchanged. The 
interaction length scale CTy between any two small A par¬ 
ticles is aAA = 1.0 and similarly for between one small A 
and one large B particles is ctas = 1.18. The interaction 
length scale between two large B particles is asB = 1-4. 

We have also done simulations with cr^ varied so that 
the system traverses from a crystalline state to an amor¬ 
phous state. This is known as amorphization transition 
[24]. Here, an effective interaction length scale Ue// 
is maintained such that aeff = , where 

xi = a ;2 = 1/2 are the fractions of the two components. 
A = < 1 measures the degree of bi-dispersity in the 

system, an and 022 are determined using the condition 
that cTe// is the same as that for original model described 
before. We kept Tco = 1.60aij. The simulations are done 
in both two and three dimensions to study the effect of 
dimensionality in the results reported. 

NVT MD simulations are done in a cubic simula¬ 
tion box with periodic boundary conditions in both two 
and three dimensions for all the model systems using 
LAMMPS package and visualization of the MD trajec¬ 
tory is done using VMD open source software [26, 27]. 
We use the modified leap-frog algorithm with the Berend- 
sen thermostat to keep the temperature constant in the 
simulation runs as implemented in LAMMPS. Any other 
thermostat does not change the results quantitatively as 
we are mostly interested in configurational changes in 
the system instead of momentum correlations. Length, 
energy an d time sca les are measured in units of aAA, 
caa and \Ja'Aj^l caa- The integration time steps used is 
dt = 0.005 in the studied temperature range. The num¬ 
ber of particles used for is iV = 2000 and the bi-dispersity 
ratio was 50 : 50. 
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FIG. 2. The schematic picture of a patchy spherical colloidal 
particle (red circle) with patches (blue and green semi-circles) 
along the equator. 


B. Model B: 

While Model A, described above, has isotropic interac¬ 
tions, in order to study whether the results reported here 
are generic for all glass forming liquids, we also perform 
simulations with a model solid with strongly anisotropic 
interactions viz. a model colloidal solid with patchy in¬ 
teractions in two dimensions. The details of the poten¬ 
tial modelling the patchy colloid is based on the angle- 
dependent interactions of Hamaker[28, 29] and Lennard- 
Jones (LJ) type as implemented in the MD simulation 
package LAMMPS[26]. In our model, each molecule con¬ 
sists of one central, large, spherical particle with six small 
equidistant patches of alternating types on it’s equator. 
In Fig.2, we draw a schematic for the patchy colloid par¬ 
ticle. We refer to the central particle shown as a large red 
circle as a T ype — 1 particle and the blue and green semi¬ 
circular patches on it’s equator as Type — 2 and Type — 3 
particles respectively. 

The interaction between two Type — 1 particles is a 
Hamaker interaction [28] between two large size colloid 
particles. The interaction between a type — 1 particle and 
a type — 2/type — 3 particle is the interaction between a 
large size colloid particle and a solvent LJ particle. Two 
type — 3 particles or a type — 2 particle and a type — 3 
particle interact with simple LJ interaction. Details of 
the interactions are given in the appendix of [30]. The 
sizes and interaction strengths between type — 1, type — 2 
and type — 3 particles are also tabulated in the appendix 
of [30]. All quantities in the table are expressed in re¬ 
duced units. The unit of length and energy are tru and 
All respectively. We choose the mass of each molecule 
ail = All = m = 1 without loss of generality, where m 
is the mass of each particle. 

We carry out molecular dynamics (MD) simulations on 
this system with an integration time-step 6t = 5x 10“^ in 
constant NAT ensemble, where N = 864 is the number of 


patchy-colloids and A is the area in two dimensions. The 
system-temperature is kept fixed at the desired value us¬ 
ing dissipative particles dynamics (DPD) thermostat [31] 
as implemented in LAMMPS. 

C. Deformation protocols: 

In all these systems the numerical experiments are de¬ 
signed as follows. We first equilibrate the model liquids 
at some high temperature and then cool it to low tem¬ 
perature below the experimental glass transition temper¬ 
ature defined as the temperature where the relaxation 
time becomes 10®. Then at that low temperature we 
perform constant pressure and temperature(NPT) simu¬ 
lations at zero pressure such that one can now remove the 
periodic boundary condition. The initial simulations at 
high temperature and high pressure was done for 5 x 10® 
time steps, then the system temperature was gradually 
reduced in 5 x 10® steps, and finally the NPT simulation 
at zero pressure was run for 5 x 10® steps. We now de¬ 
fine two side walls at the two ends of the solid in the x 
direction by pinning the particles in the end region. The 
typical size of this wall is around three inter-particle di¬ 
ameter. The other boundaries are made free. Next the 
walls are moved by an increment equal in size and op¬ 
posite in sign, i.e., the system is subjected to uniaxial 
strain. If the material is brittle it cannot deform much 
before it fractures, while for a ductile substance the dis¬ 
tance between the walls increases almost by a factor of 
two before a thin neck forms due to the plastic deforma¬ 
tions and eventually breaks. To quantify ductility of the 
system we use the percent of elongation of the system 
before breaking as given by 

r = X 100 (2) 

where, Lq is the initial (before pulling) length of the sys¬ 
tem along the tension direction ( x direction ) and L/ is 
the final length before the system breaks into two parts. 

III. RESULTS 

In this section we present results from our simulations 
following the protocols described for each of the mod¬ 
els. The bulk of our results are for the solid described 
by Model A. We show later that our main finding, i.e. 
the close relationship between the vibrational density of 
states and fracture behaviour of solids is also borne out 
in Model B. 


A. Model A: 

In Fig. 1 the snapshots of a brittle (top panel) and 
ductile material (bottom panel) under uniaxial strain. 
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FIG. 3. Top panel: Variation of ductility parameter F (see 
text for the definition) with cut-off range for a model two di¬ 
mensional glass forming liquid. F changes somewhat sharply 
as the interaction range Vco is increased and seems to reach a 
plateau once r^o becomes large enough to include the second 
neighbour shell in the pair correlation function (see text for 
details). Bottom panel: similar analysis done for the same 
model system in three dimensions. 

clearly show that brittle and ductile materials respond 
to external strain very differently. In the case of a brit¬ 
tle material, the material cannot withstand strain for a 
large deformation and hence it fails quickly while for a 
ductile material there is a neck formation which helps 
it to withstand strain for a longer period. This behav¬ 
ior can be quantified using the parameter F as described 
in the previous section. In the upper panel of Fig. 3, 
we show the dependence of F as a function of increas¬ 
ing range of interaction of the inter-particle potential, 
Tc for the two dimensional glass forming liquids. One 
can clearly see that the ductility parameter modeled by 
F is somewhat small and does not increase much until 
the interaction range starts to include the second neigh¬ 
bours. The first neighbour shell is up to a distance 1.4 in 
units of inter-particle diameter as obtained from the dip 


of pair correlation function, g{r) after the first peak. F 
then somewhat sharply increases once rc is increased be¬ 
yond the first interaction shell and then tend to saturate 
once the interaction range increases beyond the second 
neighbour distance. Similar behaviour is observed for the 
three dimensional system as shown in the lower panel of 
Fig.3 although the increase in F is somewhat smoother 
than two dimensional system. 


X 10^ 




FIG. 4. Top panel: Variation of excess vibrational density 
of states (VDOS), g{uj)lLo‘^ as a function of lu with cut-off 
range, r^o for the model A in three dimensions. Notice the 
rapid change this quantity as a function of increasing cut¬ 
off range, Vco- Bottom panel: excess VDOS obtained for the 
same model in two dimensions. The variation with increasing 
Vco is really large. 

Another thing to notice is that the variation of F with 
the cut-off range of the inter-atomic potential in two and 
three dimensions are quite large. The value of F changes 
almost 3.3 times when we vary the cut-off from 1.2 to 2.2 
for the two dimensional system, while for the three di¬ 
mensional system the change is approximately 1.6 times. 
Hence, in both the systems we could see that as the cut¬ 
off increases the amount of plastic deformation that the 
system can withstand also increases dramatically before 
failure. This can be attributed to the variation in the 
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number of neighbors of the atoms as we change the cut¬ 
off. Cut-off value of 1.2 corresponds to the first peak in 
the radial distribution function for the system. So in this 
case, the potential range is shorter than the first shell of 
the neighbors. Hence in this case, while the system is 
subjected to an external strain if any cavity forms which 
is larger than the first shell of neighbors it will not heal 
since there will be no attractive forces across it leading 
to abrupt failure. As we increase cut-off range more and 
more neighbors come in and we need a larger cavity to in¬ 
duce failure [13] or the mechanism of failure may change 
completely. 

To understand what actually changes in the system 
with increasing interaction range that leads to the ob¬ 
served brittle to ductile transition, we calculated the vi¬ 
brational density of states (VDOS) to investigate a pos¬ 
sible relationship between the mechanical properties of 
the system with that of the vibrational properties of the 
system as suggested in [13, 24, 25]. The density of states 
were calculated from the eigenvalues of the Hessian ma¬ 
trix obtained at the potential energy minimum of the sys¬ 
tem explored at that temperature. The Hessian matrix is 
defined as = d'^U/drfdr^ where U is the potential 
energy of the system which is a function of the particle 
coordinate with i indicates the particle index and a is 
different components of the space dimension (x, y or z). 
The potential energy \J is first minimized with respect 
to the particle coordinates using conjugate gradient 
methods and then used LAPACK routine to diagonalize 
the matrix which is Nd x Nd matrix, where d is the 
spatial dimension. 

The resulting VDOS g (uj) obtained as a function of 
w is reduced to the form of g (uj) where d is the 

dimensionality, and the reduced VDOS is plotted in Fig. 
4. This is done so to plot only the excess part of the den¬ 
sity of states over Debye density of state of ideal solid 
where g (uj) ~ One can clearly see in Fig. 4, 

that the excess VDOS increases dramatically with in¬ 
creasing interaction range both in two and three dimen¬ 
sional model system. This strongly suggests that excess 
vibrational density of state intimately connected the me¬ 
chanical properties of these systems as also suggested in 
Ref. [13]. 

The correlation between the excess VDOS and ductil¬ 
ity of materials although suggest interesting connection 
between the low energy vibrational states with plasticity 
in the system but does not prove the unique dependence 
of ductility on excess VDOS. To understand whether ex¬ 
cess VDOS can uniquely determine the ductility of any 
material, we have studied the brittle-ductile transition 
across the amorphization transition where the ratio of 
particle diameter of the binary glass forming liquids is 
changed systematically from 1. This way one can go from 
crystal to glassy state by tuning the ratio of particle di¬ 
ameter in the model. Notice that as we will be going 
from crystal to a glassy state with varying the diameter 
ratio, the excess VDOS will increase very sharply as an 
ideal crystal will not have any excess VDOS and it will 




FIG. 5. Top panel: Variation of F with size disorder charac¬ 
terized hy R — gaaIo-bb (see text for details) for model A 
in two dimension. Notice the sharp change in F at i? ~ 0.85. 
Bottom panel: Similar analysis done for the same model in 
three dimensions. The sharp decrease in F at i? ~ 0.90 is 
clearly visible. 


gradually increase as one goes to amorphous solids. This 
is an ideal set up to understand the role of excess density 
of state in brittle ductile transition as the excess VDOS 
can be tuned very systematically. 

In Fig. 5 we have shown the variation of F with differ¬ 
ent values ofthe ratio ofparticle diameter, R = aAA/o'BB 
for two and three dimensional systems respectively, ctaa 
is the diameter of the A type particle and asB is the 
diameter of B type particle in a typical binary mixture. 
Size ratio R = 1 corresponds to a mono-disperse system 
which will crystallize if the temperature is decreased be¬ 
low the freezing temperature. R is varied from 1.0 to 
the lowest size ratio of 0.714 is similar to the size ratio 
used in the model system mentioned above. The inter¬ 
action range is = 1.60. In the top panel of Fig.5, the 
variation of F is shown as a function of R for two dimen- 
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FIG. 6. Top left panel: Variation of VDOS as a function of uj with size disorder R for model A in two spatial dimensions. 
Notice that g{io) at small uj gains more weight as one goes from crystal to amorphous solids via the amorphization transition 
(see text for details) by decreasing R. Top right panel: g{uj)/uj is plotted to clearly show the increase in excess VDOS at small 
UJ. Bottom panels: Same analysis done for the model A in three dimensions. 


sional system and one can see that F remains more or 
less independent of i? up to i? < 0.85 and then sharply 
decreases at i? ~ 0.85 and remains independent again 
up to i? = 1.0. r changes by a factor of 7.0 in this 
model system. Similar behaviour is found for the three 
dimensional case with F changes by a factor of 2.0. It is 
important to notice that the change of F as a function of 
R in both dimensions, remains very sharp, suggesting a 
possible sharp transition from brittle to ductile materials 
at a critical dispersity ratio R = ~ 0.85. 

Now to prove the unique connection between the excess 
VDOS and the ductility we have calculated the VDOS for 


all these systems which are plotted in Fig.6. In top left 
panel of Fig.6, we have shown the VDOS for the two di¬ 
mensional model system for different dispersity ratio R 
and in right panel of the same figure the excess VDOS 
is plotted. One can clearly see that for R > 0.85, the 
excess VDOS is almost negligible as g{u!)/uj is almost in¬ 
dependent of UJ in the small uj limit and increases sharply 
once R < 0.85 and becomes almost insensitive to change 
in R for lower values. Similar behaviour is observed for 
the three dimensional case also. 

These observation clearly establishes the one to one 
connection between the excess VDOS and the ductility 
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FIG. 7. Left panel: Variation of excess VDOS, g{uj)/u} as a function of uj for the crystalline state (red circles) and the quenched 
amorphous state (blue triangles). Right top panel: Snapshot picture of the material in crystalline state under uniaxial strain 
along x-direction. Similar snapshot in the amorphous state (bottom right panel). 


of a material irrespective of whether it is a crystalline 
solid or an amorphous solids. As mentioned earlier, in 
[25], it was shown that mechanical properties of a slightly 
defected crystal with a very small density of vacancies is 
more close to a jammed amorphous solid than to the 
perfect crystal. They also pointed out that vibrational 
density of states of a solid with even a very high degree of 
crystallinity resembles more closely to that of the amor¬ 
phous jammed solid. A sharp transition in mechanical 
and vibrational properties is also suggested in that work. 
Our work here similarly suggest a possible sharp brit¬ 
tle to ductile transition as one goes from crystalline to 
amorphous solid by changing dispersity parameter, R. It 
further points out a one-to-one correspondence between 
the mechanical and vibrational properties of solids. 

A very similar connection between the observed be¬ 
haviour in our work and behaviour observed in fiber bun¬ 
dle model [32, 33] with increasing dispersity in the fiber 
strength is worth mentioning. The fiber bundle model is 
a simplistic model to understand fracture in materials. 
In this model, two plates are connected by fiber bundles 
and then the plates are loaded to mimic the tension ex¬ 
periments [32, 33]. The solid will then be modeled via 
the fiber in between the plates. Now if one takes the 
strength of these fibers to be same then it will model a 
homogeneous solid, but if one takes the strength of the 
fiber to be heterogeneous, then it will be closely mimick¬ 
ing the behaviour of an amorphous solids. In [34], one 
such model is studied with the strength of fibers taken 
from a Gaussian distribution of width 5. This disorder 
parameter (5, was then varied and the mechanical prop¬ 
erties of the solid is studied. It was shown that as one 
increases the disorder parameter 5, the solid modeled by 


the fiber bundle model shows a brittle to ductile transi¬ 
tion at a critical disorder strength S = Sc- Our results 
seems to strongly suggest such a transition in realistic 
model materials [35], thereby providing a nice model sys¬ 
tem to test the prediction of these hber bundle models for 
further improvement of these minimalistic models [36]. 


B. Model B: 

So far we have tested the connection between the ex¬ 
cess VDOS and the ductility for system where the inter¬ 
action potential between two particles is isotropic. To 
test whether the same connection still holds for other 
model interaction potentials with anisotropic interaction, 
we have performed the similar analysis for a patchy col¬ 
loid model described in the model and simulation de¬ 
tails section before. This model at different temperature 
and density shows many interesting phases (see [30] for 
further details and the phase diagram), e.g triangular, 
honeycomb, square lattice and their corresponding phase 
coexistence. For a large portion of the phase diagram, 
this model also shows amorphous structure. In this case, 
an amorphous solid may be formed by quenching rela¬ 
tively rapidly into a region of the phase diagram where 
there are many competing crystalline states. The solid, 
instead of selecting between these almost degenerate free 
energy minima corresponding to very different structures, 
chooses to remain amorphous or forms microscopic crys¬ 
tallites each with wildly varying local coordination. 

We perform loading experiments using the protocol de¬ 
scribed in section II for this model solid at two represen¬ 
tative states. Both the solids are at a number density 





























p = 0.77, where the system at low temperature forms a 
honeycomb lattice. The first state is crystalline with two 
dimensional honeycomb order obtained by arranging par¬ 
ticles in an initial crystalline order and equilibrating the 
structure at the chosen density and temperature. The 
second state, at the same density, is amorphous and is 
formed by a temperature quench from the high tempera¬ 
ture liquid. Note that this amorphous structure is clearly 
in a state of dynamic arrest. 

For the crystalline state one has no excess density of 
state over the Debye theory and for the amorphous state 
there is an excess of states as shown in Fig.7a. Now, ac¬ 
cording to our previous observations, the state with ex¬ 
cess VDOS should be more ductile than the one with less 
excess VDOS. It is indeed true that the amorphous state 
with excess VDOS can withstand longer tensile strain 
and bears F ~ 32, while the pure crystalline state with 
less VDOS breaks like a brittle materials with F ~ 1.5. 
Fig.7b shows one instance of fracture each in the crys¬ 
talline and amorphous states obtained for the patchy col¬ 
loidal system with different preparing protocol. The crys¬ 
talline state shows a sharp fracture while the amorphous 
state seems to show more heterogeneous fracture profile. 
The increase in low frequency vibrational states for the 
amorphous structure is quite dramatic. Thus we believe 
that the excess VDOS uniquely captures the mechanical 
behaviour of a solid with any form of inter-atomic in¬ 
teraction potential and irrespective of whether it is in a 
crystalline or amorphous state. 


IV. CONCLUSIONS 

In conclusion, we have done extensive MD simulation 
of different glass forming liquids with both isotropic and 
anisotropic pairwise interactions to understand the mi¬ 
croscopic origin of ductility in amorphous solids. There 
has been many attempts to link macroscopic mechani¬ 
cal properties of solids such failure mechanism with mi¬ 
croscopic “atomic” interactions. For example, extensive 
calculations and analysis of known experimental data for 
a number of metallic and non-metallic solids by Rice 
and Thomson [37] shows that ductility in solids is re¬ 
lated to the ability of a crack tip to nucleate dislocations. 
Such dislocations nucleated at crack tips blunt the crack 


and produce ductile behaviour. A confirmation of this 
scenario has been found recently in the computer sim¬ 
ulations [38, 39] where dislocation tangles surrounding 
cracks in ductile materials is readily observed in large 
scale computer simulations. 

In this paper, on the other hand we find a different 
kind of correspondence. We observe that ductility of a 
material is intrinsically connected to excess density of 
vibrational states over the Debye density of state of an 
ideal solid. Findings of Goodrich et.al. [25], clearly sug¬ 
gest that excess vibrational density of state is a unique 
characteristic of the material which strongly connects to 
the mechanical property of the system. They also showed 
that solid with not perfect but large crystalline order 
show properties closely matching with that of an amor¬ 
phous jammed solid than ideal crystal with defects. The 
vibrational density of state of these systems becomes very 
close to that of an amorphous solids. Here we showed 
that excess vibrational density of states uniquely deter¬ 
mines the ductility of an amorphous solids under exter¬ 
nal tensile load. We also showed that there exists a sharp 
brittle to ductile transition as one drives the system from 
crystalline to amorphous state via amorphization transi¬ 
tion which is in close agreement with the predictions of 
minimalistic fiber bundle model [34]. It will be nice in 
future to understand in details this transition by doing 
extensive finite size scaling analysis. Study along these 
directions are in progress and will be published else where 
[36] Our results on the patchy colloidal model also con¬ 
firms that the reported results are rather generic for any 
glass forming liquids. 

Is there a connection between VDOS and the ability 
of the solid to produce dislocations? While a complete 
answer to this very interesting question has not been 
found, there are intriguing possibilities. In [40] it was 
shown that probability of formation of defects in crys¬ 
talline solids depends on the amplitude of certain non- 
affine displacement fluctuations [41] which act as precur¬ 
sors. Non-affine defect precursors also cause instabilities 
in the crystal leading to mode softening which may, in 
turn enhance the VDOS in the small frequency regime. 
We thus feel that our finding of direct connection be¬ 
tween excess density of vibrational states and plasticity 
in amorphous solid as well as crystalline solids will help 
us better understand plasticity in these materials within 
an unified theoretical framework. 
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